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Abstract 

Gene exchange between species occurs in areas of secondary contact, where two 
species have the opportunity to hybridize. If heterospecific males are more 
common than conspecific males, females will experience more encounters with 
males of other species. These encounters might increase the likelihood of het- 
erospecific matings, and lead to the production of hybrid progeny. I studied 
the mating behavior of two pairs of sibling species endemic to Africa: Drosoph- 
ila yakuba/Drosophila santomea and Drosophila simulansi Drosophila sechellia. 
Drosophila yakuba and D. simulans are cosmopolitan species widely distributed 
in the African continent, while D. santomea and D. sechellia are island endem- 
ics. These pairs of species hybridize in nature and have the potential to 
exchange genes in natural conditions. I used these two pairs of Drosophila spe- 
cies, and constructed mating communities of different size and different hetero- 
specific:conspecific composition. I found that both the total number of 
potential mates and the relative frequency of conspecific versus heterospecific 
males affect female mating decisions in the cosmopolitan species but not in the 
island endemics. These results suggest that the population characteristics, in 
which mating occurs, may affect the magnitude of premating isolation. Com- 
munity composition might thus facilitate, or impair, gene flow between species. 



Introduction 

When different animal species come together in the 
same geographical location and share at least part of 
their habitat, biological traits associated with mate 
choice can prevent interbreeding (Coyne and Orr 2004; 
Price 2007). One of the most effective mechanisms of 
reproductive isolation results from mate preferences 
that differ between species (Mayr 1942; Kirkpatrick and 
Ravigne 2002; Coyne and Orr 2004; Ritchie 2007). 
Premating behavioral isolation occurs when one or both 
partners discriminate against the other, thus precluding 
mating and gene flow (Kaneshiro 1980; Safran et al. 
2013). Because of their greater investments of resources 
during and after mating, females usually are the ones 
that exert the choice and discriminate against hetero- 
specific males. Mate choice can vary in response to fac- 
tors intrinsic to the potential partners (i.e., age, health, 
reproductive fitness), or to extrinsic factors in the envi- 
ronment in which potential mates encounter each other 
(i.e., habitat and phenology differences; Rice and Hostert 



1993; Rolan-Alvarez and CabaUero 2000; Knowles et al. 
2004; Coyne et al. 2005). 

When heterospecific matings occur, and there is not 
complete intrinsic reproductive isolation between the spe- 
cies, interspecific hybrids are produced. The consequences 
of hybridization can vary. One potential outcome is that 
the production of unfit hybrids might impose selection to 
make premating isolation stronger. In this process, 
termed reinforcement, the enhancement of premating iso- 
lation occurs as a byproduct of selection against maladap- 
tive hybridization (reviewed in Servedio and Noor 2003). 
A second outcome is that gene flow between species 
might lead to introgression of foreign genes into the gen- 
ome; gene flow can have deleterious consequences (e.g.. 
Fang et al. 2012), but also may become raw material for 
the origin of new adaptations (reviewed in Hedrick 
2013). Finally, in rare instances, the production of hybrids 
can lead to the origin of a new lineage that shows repro- 
ductive isolation toward the parentals and thus consti- 
tutes a new hybrid species (Arnold and references therein, 
Schumer et al. 2014). 
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Even though there is consensus about the evolutionary 
importance of hybridization in nature, there is not a 
comprehensive understanding of what biotic or abiotic 
factors might facilitate hybridization between species. Sev- 
eral arguments posit that if heterospecific males are much 
more common than conspecific males, then females will 
experience more encounters with heterospecifics, and pro- 
duce more hybrid progeny than females in areas where 
conspecifics are common and heterospecifics are rare 
(Volpe 1955; Waage 1979; Peterson et al. 2005; Nosil 
2012). Nonetheless, the assumption that differential rela- 
tive frequency of conspecific and heterospecific males can 
lead to different levels of hybridization has remained 
untested. Here, I used two pairs of Drosophila species to 
assess whether different ratios of conspecifics to hetero- 
specifics can lead to increased hybridization. 

Drosophila yakuba and D. santomea are two sister spe- 
cies of the melanogaster subgroup of species (Lachaise 
et al. 2000). Drosophila yakuba is a human-commensal 
species that is widespread throughout sub-Saharan Africa 
(Burla 1954, Lachaise et al. 1988). Drosophila santomea, 
the sister species of D. yakuba, is an endemic species to 
the highlands of Sao Tome, a volcanic island off the coast 
of Cameroon (Lachaise et al. 2000). On Pico de Sao 
Tome, D. yakuba inhabits low elevations (below 1450 m), 
and is found in open habitats (Llopart et al. 2005a). In 
contrast, D. santomea inhabits the mist forests of the 
island at elevations between 1153 m and 1800 m (Llopart 
et al. 2005a, Llopart et al. 2005b). Even though these two 
species diverged around million years ago, they occasion- 
ally hybridize in the midlands of the mountain Pico de 
Sao Tome in a recent area of secondary contact (Llopart 
et al. 2005a, Llopart et al. 2005b). Hybrids between these 
two species show reduced fitness as a result of intrinsic 
postzygotic isolation: hybrids of both sexes show mildly 
reduced viability when compared to pure species. 
Hybrid females are usually fertile but Fi hybrid males, 
and a large proportion of the males from advanced inter- 
crosses, are sterile and show gross defects in spermatogen- 
esis (Moehring et al. 2006, Matute et al. 2009; Matute 
and Coyne 2010). 

Drosophila simulans is a human commensal that is 
thought to have originated in Madagascar and currently 
has a cosmopolitan geographic range, including the 
Seychelles archipelago in the Indian Ocean. Drosophila se- 
chellia (Tsacas and Biichli 1981; David et al. 2007), an 
endemic species to the Seychelles, is closely related to 
D. simulans (Amlou et al. 1998; Kliman et al. 2000; 
McDermott and Kliman 2008). The two species diverged 
approximately 0.2 million years ago (Kliman et al. 2000; 
McDermott and Kliman 2008). When the two species are 
crossed, they produce fertile hybrid females but sterile -Fi 
males (Coyne and Kreitman 1986; Coyne et al. 1991; 



Cabot et al. 1994; HoUocher and Wu 1996). In a similar 
manner to the D. yakuba/D. santomea species pair, the 
majority of males from advanced intercrosses are sterile 
as well (Moehring et al. 2005). Drosophila simulans and 
D. sechellia currently hybridize in the central islands of 
the Seychelles archipelago (Matute and Ayroles 2014) and 
are thought to have experienced gene flow at some point 
(Solignac and Monnerot 1986; Garrigan et al 2012), indi- 
cating that hybridization has played a role in their evolu- 
tionary history. 

In both species pairs, and both directions of the cross, 
heterospecific matings occur at a much lower frequency 
than conspecific matings. In nonchoice experiments (i.e., 
one female and one male in a vial), females from both 
D. yakuba/D. santomea and D. simulans/D. sechellia spe- 
cies pairs show a strong tendency to mate with conspe- 
cific males and rarely accept heterospecific males (Coyne 
1992; Coyne et al. 2002). 

In Drosophila, attributes of individuals and their envi- 
ronments may affect the frequency of hybridization but 
they remain largely unstudied in experimental settings. 
Coyne et al. (2005) explored the effects of three factors - 
individual female health, the effect of female having more 
space to escape from unwanted mates, and the effect of 
having a fruit present in the mating arena - and found 
that none of them had an effect on female mating deci- 
sions. Many other attributes of the mating community 
remain to be tested. 

Population composition may affect female mate choice 
and the frequency at which hybridization occurs. Here I 
show that the two population factors, the relative ratio of 
heterospecific versus conspecific males, and the total 
number of potential mates, affect female mating decisions 
in the aforementioned Drosophila species pairs. Gene flow 
between species may therefore be contingent upon the 
communities in which members of those species encoun- 
ter one another. 

Methods 

Stocks 

All collected stocks and populations were reared on stan- 
dard cornmeal/Karo/agar medium at 24°C under a 12-h 
light/dark cycle. I used genetically heterogeneous strains 
of each species (i.e., synthetic lines) by combining virgin 
males and females from several isofemale lines collected 
in Sao Tome outside the hybrid zone (i.e., aUopatric 
lines). For all the experiments involving D. yakuba and 
D. santomea, I used the D. yakuba SYN2005 and 
D. santomea SYN2005 stocks, respectively (Matute et al. 
2009; Matute and Coyne 2010). For the experiments 
involving D. simulans and D. sechellia, I used D. simulans 
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Florida city (Coyne and Beecharn 1987) and D. sechellia 
SYNDenis (Matute and Ayroles 2014) synthetic stocks. All 
stocks were kept in large numbers after they were created. 

Effect of relative frequency on mating 
behavior 

I explored two components of mating behavior, conspe- 
cific copulation latency (i.e., how long does it take for a 
conspecific mating to take place) and conspecific copula- 
tion duration (i.e., how long did copulation last) when 
Drosophila females of four different species (D. yakuba, 
D. santomea, D. simulans, and D. sechellia) were exposed 
to mating situations in which conspecific males were 
present and males from a second hybridizing species were 
also present. I assayed all the possible combinations 
between the two factors: the relative frequency of hetero- 
specific males (10 different frequencies from 10% to 99% 
in increments of 10%) and the total number of flies in 
the mating assay (10 different totals from 100 to 1000 in 
increments of 100). In total, I recorded both premating 
isolation behaviors in 100 combinations and assayed 15 
females per combination. 

All measurements of premating isolation were carried 
out using previously described methods (Coyne et al. 
2002; Matute 2010); unlike previous experiments, females 
were allowed to choose their mate. Briefly, females and 
males were collected as virgins within 8 h of eclosion. All 
individuals were kept in vials of 22 flies of the same sex for 
3 days. On Day 3, females were housed in red-colored 
food (which turned their abdomens red) for their easier 
identification in the experimental mating population. Red 
dye has a negligible effect on mating choice (Ting et al. 
2001; Matute 2013). Four days after collection, a single vir- 
gin female and the virgin males dictated by each experi- 
mental combination, as described above, were combined. 
All flies were transferred without anesthesia to a bottle 
with cornmeal food. For each mated female, I recorded 
conspecific copulation latency (i.e., how long did it take 
for a D. yakuba female to find a conspecific mating partner 
and start copulating), and conspecific copulation duration 
(i.e., how long did conspecific females and males remained 
joint in mating). Flies were observed for 1 h. I did not 
observe any heterospecific matings using this approach. 

In parallel, I set up mating trails in which the only 
males present were conspecific. I varied the number of 
flies in the mating assay (10 different totals from 100 to 
1000 in increments of 100). I recorded both conspecific 
copulation latency and conspecific copulation duration in 
mating trials of ten different sizes and assayed 15 females 
per population size. 

To study heterospecific mating frequency (i.e., how 
often females accepted heterospecific males even though 



conspecific males were present in the population), I mea- 
sured premating isolation over 24 h for a single female 
that were housed with both conspecific and heterospecific 
males in the numbers and relative frequencies described 
above. I let trials proceed for 24 h and then anesthetized 
all the flies in the mating population with CO2, discarded 
the males and returned each female to the vial where the 
mating took place in order to observe the resulting Fi 
generation. I started 100 replicates per combination of 
variables (i.e., each combination of heterospecific relative 
frequency and total number of flies). I raised the progeny 
of each female by standard fly husbandry methods, and 
used the presence or absence of hybrid progeny as a con- 
servative proxy of whether heterospecific matings 
occurred. The frequency of heterospecific matings for 
each treatment was calculated by counting how many of 
the vials within a block produced hybrid progeny. For 
D. yakuba and D. santomea experimental populations, I 
used abdominal pigmentation to identify whether vials 
produced hybrid progeny. I qualitatively scored the 
abdominal pigmentation of 100 males per vial. Drosophila 
yakuba males have a dark abdomen, while D. santomea 
males have a light abdomen. _Fi males have an intermedi- 
ate abdominal pigmentation (Llopart et al. 2002). For 
D. simulans and D. sechellia, experimental populations, I 
used male genital morphology to identify whether vials 
produced hybrid progeny. I qualitatively scored the mor- 
phology of the posterior genital arch of 100 males per 
vial. Drosophila simulans males have red spheroid large 
genital arches. Drosophila sechellia males have small and 
elongated genital arches. The hybrids show intermediate 
genital morphology and can be distinguished from the 
parental species (Coyne and Kreitman 1986; Coyne et al. 
1991; MacDonald and Goldstein 1999; Matute and 
Ayroles 2014). 

Statistical analyses 

AU statistical analyses were carried out using R (R Devel- 
opment Core Team 2005). 

To determine the significance of the frequency of het- 
erospecific males and the mating population size on mat- 
ing behavior, I analyzed the conspecific copulation 
latency and conspecific copulation duration, I fitted a 
multiple regression for each of the two components of 
premating behavior. The premating behavior trait was the 
response, while the relative frequency of heterospecific 
males and the total number of flies were the continuous 
variables. Multiple regressions took the form: 

Yij = Fi + S, + Fi xSj+ E,j 

where Yy is the response (behavioral trait), f,- is the fre- 
quency of heterospecific males, Sj is the mating population 
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size (i.e., number of males in the mating assay), f,- x Sj is 
the interaction between the variables, and _Ey is the error 
term. The significance of independent relationships in the 
multiple regressions was determined using a two-tailed t- 
test on the regression coefficients (degrees of freedom, 
df = 1496). 

For the heterospecific mating frequency data (collected 
when females were housed with different combinations of 
males over 24 h), I used the "aod" package in R (Lesnoff 
and Lancelot 2012). I fitted generalized logistic models 
with binomial response distributions in which whether a 
vial produced hybrid progeny or not was the response. 
The total number of flies in the vial and the relative fre- 
quency of heterospecifics were considered continuous 
variables. I allowed for an interaction term between these 
two variables. Significance of each variable was assessed 
using a Wald text following a distribution (df =1). 
Linear models followed the form: 

Yij = Fi + Sj + F, X Sj + E,j 

P values under 0.05 were considered significant. 

I assessed whether total number of flies and the ratio 
of heterospecifics were multicoUinear variables (i.e., the 
two predictor variables were correlated). The reason for 
this concern was that the ratio of heterospecifics was cal- 
culated using the total number of flies as the denomina- 
tor. This can lead to autocorrelation of the two variables. 
I calculated the variance inflation factor (VIF) between 
the two regression coefficients using the R package HH 
(Heiberger 2009). The VIF for predictor i equals: 



where is the coefficient of determination from a 
regression of predictor i against the remaining predictors. 
In the case of the eight multiple regressions shown in this 
study, we calculated the VIF values for the two predictor 
variables in models that had no interaction term. Because 
all the multiple regressions had the same two predictor 
variables (and the same values), I only needed to calculate 
one VIF. The VIF for the model 

R,j = Fi + Sj + Eij 

where Rjj is either latency or duration, Fj is the ratio of 
heterospecifics and is the mating population size, 
equaled 1. Values of VIF larger than 5 are considered evi- 
dence of coUinearity between predictor variables (Stine 
1995; O'Brien 2007). This result indicates that the two 
predictor variables are not multicoUinear. 

To plot heterospecific mating frequency data, I used 
Akima interpolation splines (Akima 1970) between sam- 
pling intervals, and plotted them using contour maps 
with the "akima" R package (Akima et al. 2006). 



D. R. Matute 

Results 

D. yakubalD. santomea 

I studied the mating sexual behavior of D. yakuba females 
exposed to different relative frequencies of D. santomea 
and D. yakuba males in populations with different total 
numbers of males. The mating behaviors I studied were 
conspecific copulation latency (time until mating began) 
and duration. The average conspecific copulation latency 
and conspecific copulation duration per treatment are 
shown in Figure 1. Both the relative frequency of 
D. santomea males and population size affected the mat- 
ing behavior of D. yakuba females toward conspecific 
males; in larger populations, the conspecific copulation 
latency was longer (Table 1, Fig. 1) and conspecific copu- 
lation duration was shorter (Table 2, Fig. 1). In the case 
of duration, but not of latency, the interaction between 
population size and composition was also significant, 
indicating a strong interplay between the size of the mat- 
ing community and the ratio of heterospecific to conspe- 
cific males (Tables 1, 2). 

Notably when a single D. yakuba female was housed 
with only conspecific males in populations of different 
sizes, both copulation latency and duration remained 
constant (Table 3). These results indicate that a simple 
increase in the population size (and consequently the sex 
ratio) does not lead to changes in mating behavior. Such 
changes are only observed when there are heterospecific 
males present in the assay. 

In addition to conspecific copulation latency and dura- 
tion, I investigated the effect of community composition 
on reproductive isolation in D. yakuba. There were no 
heterospecific matings observed in the 1-h trials, so I next 
looked at heterospecific mating frequencies in 24-h trials 
when D. yakuba females were housed with the same ratios 
of conspecific:heterospecific males and total numbers of 
males as in the 1-h observations. In this experiment, large 
populations that had a high ratio of heterospecific to con- 
specific males were more prone to produce hybrid prog- 
eny (binomial logistic regression, Wald test, = 384.9, 
df = 3, P < 1 X 10"'°, Fig. 2). I found that the ratio of 
heterospecific to conspecific males had a significant effect 
on the frequency of heterospecific matings (Wald test, 

= 59.1, df = 1, _P < 1.1 X 10"'°). The total number of 
flies also had a significant effect on the frequency of het- 
erospecific matings {-/^ = 35.4, df = 1, P = 2.9 x 10"'). 
The interaction between these two variables was also sig- 
nificant (x^ = 11.9, df = 1, P = 5.60 X 10"*). In general, 
heterospecific matings were observed only in populations 
in which more than 70% of the males were D. santomea. 

I next looked at the mating behavior of females from 
the other species of this pair. Using the same procedures 
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(A) Copulation latency 



(B) Copulation duration 



Figure 1. High ratios of Drosophila santomea 
males relative to Drosophila yakuba males 
affect two components of mating sexual 
behavior (conspecific copulation latency -CCL- 
and conspecific copulation duration -CCD-) in 
D. yakuba females. (A) The number of 
D. santomea males (%) significantly affected 
conspecific copulation latency in D. yakuba 
females. (B) The number of D. santomea males 
(%) significantly affected conspecific 
copulation duration in D. yakuba females. (C) 
The size of the mating population significantly 
affected conspecific copulation latency in 
D. yakuba females. (D) The size of the mating 
population significantly affected conspecific 
copulation duration in D. yakuba females. Each 
circle represents the average of 15 replicates. 
All significance values are shown in Table 1. 
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but with the species reversed, I investigated whether con- 
specific mating behavior or levels of reproductive isola- 
tion in D. santomea females were also affected by the 
ratio of heterospecific versus conspecific males and the 
size of the mating population. Neither of the two factors 
that I manipulated (as above, ratio of heterospecific to 
conspecific males and the mating population size) led to 
changes of the magnitude of reproductive isolation in 
D. santomea females toward D. yakuba males, nor to 
changes in the conspecific copulation latency (Table 1) or 
duration of D. santomea females (Table 2). In controls 
that only included different numbers of D. santomea 
males I observed no differences in copulation latency or 
copulation duration among treatments (Fig. 3, Table 3). 
These results indicate that the mating behavior of 
D. santomea females is not affected by the presence of 
heterospecific males regardless of the frequency at which 
they are present. Furthermore, in populations where both 
D. yakuba and D. santomea males were present, I 
observed no heterospecific copulations, neither in the 1-h 
measurements (« = 1500 observations and 100 different 
treatments) nor in the 24-h trials (n = 10,000 observa- 
tions and 100 different treatments). Although they will 
occasionally mate with D. yakuba males in no-choice 
matings (Coyne et al. 2002), D. santomea females will 
always choose D. santomea males when given the chance. 

D. simulanslD. sechellia 

In order to assess whether the results observed for D. yak- 
uba/D. santomea were specific to that species pair or 
whether they also applied to other species; I followed the 



same experimental procedures in a second species pair: 
D. simulans and D. sechellia. I measured the magnitude of 
the two conspecific mating behavior traits in D. simulans 
females exposed to different ratios of D. sechellia to 
D. simulans males. For 1-h experiments, multiple regres- 
sions showed that the frequency of the heterospecific 
D. sechellia, and the total number of flies affected both 
components of mating behavior in each mating trial for 
D. simulans females (see Methods, Table 1 and 2, Fig. 4). 
The interaction between these two variables was margin- 
ally significant for conspecific latency (Table 1) and con- 
specific duration (Table 2). I also assayed mating 
behavior in cages that contained only D. simulans males 
and observed no variation in latency or duration 
(Table 3). These results, similar to the ones from D. yak- 
uba, suggest that a simple increase in population size is 
not enough to modify the mating behavior of D. simulans 
females and that changes in mating behavior only occur 
when there are heterospecific males present in the assay 
(Table 3, Fig. 4). 

While both latency and duration of conspecific matings 
were affected by the presence and number of heterospecif- 
ics, no heterospecific matings were observed in these 1-h 
trials either. Next I examined heterospecific matings when 
I housed D. simulans females with different combinations 
of D. simulans and D. sechellia males for 24 h. Similar to 
the observations of D. yakuba, large populations with few 
conspecific males were more likely to produce hybrid 
progeny (binomial logistic regression, Wald test, = 
361.1, df = 3, P < 1.4 X 10"'^ Fig. 5). I found there was 
a strong effect on the number of heterospecific matings of 
the relative frequency of conspecific to heterospecific 
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Table 3. Mating assays witli only conspecific males and different population sizes. 



Species 


Copulation latency 






Copulation duration 






Coefficient 


t value 


P-value 


Coefficient 


t value 


P-value 


Drosophila yakuba 


-8.009 X 10"* 


-0.537 


0.592 


3.25 X 10"* 


0.144 


0.885 


Drosophila santomea 


4.746 X 10"* 


0.381 


0.704 


-4.875 X 10"* 


-0.318 


0.751 


Drosophila simulans 


2.142 X 10"^ 


1.30 


0.196 


-1.289 X 10"^ 


-0.734 


0.464 


Drosophila secheiila 


9.093 X 10"* 


0.758 


0.449 


-1.776 X 10"^ 


-1.201 


0.232 



No heterogeneity was detected in any of the four species at the two mating behaviors. I found no significant values (P < 0.05, df = 1496). 



Proportion of hybridizing females, D. yakuba Prop 




Population size 

Figure 2. The ratio of Drosophila santomea males relative to that of 
Drosophila yakuba males affects the strength of premating isolation in 
D. yakuba females. In conditions where D. santomea males are 
oven/vhelmingly more common than D. yakuba males, D. yakuba 
females accept D. santomea males. Light gray indicates conditions 
where heterospecific matings were obsen/ed. 



males (x^ = 54.7, df = 1, P = 1.9 x lO"*"), the total num- 
ber of flies in the cage (y^ = 22.7, df = 1, 
P = 1.9 X 10^''), and the interaction between these two 
variables = 6.20, df = 1, P = 0.013). 

Finally, I studied the mating behavior of females 
from the second species of this pair, D. sechellia, in 
communities of different sizes and compositions. Nei- 
ther of the two experimental variables (population size 
nor relative frequency of heterospecifics) affected the 
mating behavior (conspecific mating latency or dura- 
tion) of D. sechellia females (Tables 1, 2, Fig. 6). I 
observed no changes in copulation latency or copula- 
tion duration when the assays only included D. sechellia 
males (Table 3). 

No heterospecific matings were observed with D. sechel- 
lia for either 1-h or 24-h trials. These results indicate that, 
as was the case for D. santomea, the mating behavior of 
D. sechellia females is not affected by the presence of het- 
erospecific males regardless of the ratio of heterospecific 
to conspecific males and the mating population size. 



(A) Copulation latency 



(B) Copulation duration 



Figure 3. The ratio of Drosophila yakuba 
males relative to Drosophila santomea males 
had no effect on two components of mating 
sexual behavior (conspecific copulation latency 
-CCL- and conspecific copulation duration - 
CCD-) in D. santomea females. (A) The 
number of D. yakuba males {%) did not affect 
conspecific copulation latency in D. santomea 
females. (B) The number of D. yakuba males 
(%) did not affect conspecific copulation 
duration in Drosophila santomea females. 
(C) The size of the mating population did not 
affect conspecific copulation latency in 
D. santomea females. (D) The size of the 
mating population did not affect conspecific 
copulation duration in D. santomea females. 
Each circle represents the average of 15 
replicates. All significance values are shown in 
Table 1. 
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(A) Copulation latency 



(B) Copulation duration 
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Figure 4. High ratios of Drosophila sechellia 
males relative to Drosophila simulans males 
affect two components of mating sexual 
behavior (conspecific copulation latency -CCL- 
and conspecific copulation duration -CCD-) in 
D. simulans females. (A) The number of 
D. sechellia males (%) significantly affected 
conspecific copulation latency in D. simulans 
females. (B) The number of D. sechellia males 
(%) significantly affected conspecific 
copulation duration in D. simulans females. 
(C) The size of the mating population 
significantly affected conspecific copulation 
latency in D. simulans females. (D) The size of 
the mating population significantly affected 
conspecific copulation duration in D. simulans 
females. Each circle represents the average of 
15 replicates. All statistics and significance 
values are shown in Table 1. 



Proportion of hybridizing females, D. simulans Prop 




100 200 300 400 500 600 700 800 900 1000 

Population size 

Figure 5. The ratio of Drosophila sechellia males relative to that of 
Drosophila simulans males affects the strength of premating isolation 
in D. simulans females. In conditions where Drosophila sechellia males 
are overwhelmingly more common than D. simulans males, 
D. simulans females accept D. sechellia males. 

Discussion 

This study explores how relative frequency and popula- 
tion size, individually and interacting, influence mating 
behavior in at least two species of Drosophila. The results 
shown here have two implications for our understanding 
of how premating isolation is involved in the persistence 
of potentially hybridizing species. First, they confirm that, 
in some animal species, hybridization is more likely if 
heterospecific males are disproportionally more abundant 
than conspecifics males. The relative frequency of males 



has been proposed as a factor influencing the magnitude 
of reproductive isolation in natural populations (Harper 
and Paulson 1994; Berglund 1995; Sprenger et al. 2011; 
WiUis et al. 2011, 2012; Verzijden et al. 2012). Under this 
regime, females will accept males from other species if the 
environment in which mating takes place is not condu- 
cive to finding a conspecific partner - for instance, if 
females cannot detect the few conspecific males in the 
community, or if they encounter many heterospecific 
males before a conspecific (Wilson and Hedrick 1982; 
Willis et al. 2012). These results are of particular impor- 
tance for the study of reinforcing selection, the enhance- 
ment of prezygotic isolation as a byproduct of 
maladaptive hybridization (reviewed in Servedio and 
Noor 2003; Coyne and Orr 2004; Pfennig and Pfennig 
2009). It is commonly argued that if reinforcing selection 
affects only one species in a hybrid zone, it should be the 
rarest one (Coyne and Orr 2004; Price 2007; Nosil 2012). 
The basic premise of this assertion is that the less abun- 
dant species will be more at risk of hybridizing because 
they will face more encounters with heterospecific mates. 
This claim makes intuitive sense but has remained largely 
untested. My results provide strong evidence that for 
some Drosophila species, the less abundant species in a 
community is more at risk of mating heterospecificaUy. 

Not only are D. yakuba and D. simulans less abundant 
in their respective hybrid zones, but they are more wide- 
spread outside of the hybrid zone. Drosophila santomea 
and D. sechellia, however, are endemic to their island 
habitats. Behavioral isolation is much stronger when the 
matings involve females from the island endemics than 
when they involve females from the widespread species, 
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(A) Copulation latency 



(B) Copulation duration 



Figure 6. The ratio of Drosophila simulans 
males relative to Drosophila secheilia males 
had no effect on two components of 
premating sexual behavior (conspecific 
copulation latency -CCL- and conspecific 
copulation duration -CCD-) in D. secheilia 
females. (A) The number of D. simulans males 
(%) did not affect conspecific copulation 
latency in D. secheilia females. (B) The number 
of D. simulans males (%) did not affect 
conspecific copulation duration in D. secheilia 
females. (C) The size of the mating population 
did not affect conspecific copulation latency in 
D. secheilia females. (D) The size of the mating 
population did not affect conspecific 
copulation duration in D. secheilia females. 
Each circle represents the average of 15 
replicates. All statistics and significance values 
are shown in Table 1 . 
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regardless of the experimental design used to quantify the 
reproductive isolation (Coyne 1992; Coyne et al. 2002). 
The results from this report are in the same vein of previ- 
ous studies and show that even in environments in which 
heterospecific males are overwhelmingly more common 
than conspecific males, D. santomea and D. secheilia 
females strongly discriminate against heterospecific males 
(Coyne 1996; Tomaru et al. 2004). 

Finally, I also demonstrate that the size of the mating 
population has an effect on the likelihood of the occur- 
rence of heterospecific matings. Because each assay was 
composed with a single female, when the population size 
of the male conspecific varies, the sex ratio also varies. 
This means I cannot fully disentangle the effects of popu- 
lation size and sex ratio with the current set of data, and 
direct observations of multifemale populations would 
pose distinct challenges. The results from the mating 
assays with only conspecific males do allow me to rule 
out that a simple increase in the population size (and 
consequently the sex ratio) causes changes in copulation 
latency and copulation duration. Instead, changes in these 
behaviors are only observed when heterospecific males are 
present in the assay. Regardless, these results indicate that 
the characteristics of the community in which mating 
takes place can affect the magnitude of reproductive isola- 
tion between potentially hybridizing species. 

The effects of biotic factors on hybridization are largely 
unexplored in animals but not in plants. PoUen dispersal, a 
major factor leading to gene flow in plants, is heavily influ- 
enced by heterospecific relative frequency in nature (Camp- 
bell and Halama 1993; Bosch and Waser 1999, 2001). 
Studies on pollinator competition have revealed that the 



presence of heterospecifics at different densities can affect 
the relative rates of interspecific pollen transfer (Price and 
Waser 1982; Kohn and Waser 1985; Campbell and Waser 
1989; Campbell and Halama 1993; Bosch and Waser 2001; 
reviewed in Mitchell et al. 2009). The nature of reproduc- 
tive isolation differs dramatically between plants and ani- 
mals with internal fertilization; while the results in plants 
are poUinator-mediated and thus extrinsic, the results 
shown here are intrinsic to the decision-making process of 
females choosing whether to accept a potential mate. 

Even though there is little empirical evidence for 
potential effects of the relative frequency of heterospecific 
to conspecific males on the magnitude of reproductive 
isolation in animals, there is no reason why heterospecific 
males cannot be seen as low-quality males and conspecif- 
ics as high-quality males (Penteriani 2003; Kokko and 
Rankin 2006). This is bound to be especially true in those 
cases in which postzygotic isolation is already in place, as 
happens in the two studied species pairs (i.e., hybrid 
males from the crosses D. yakuba x D. santomea, and 
from the crosses D. simulans x D. secheilia are com- 
pletely sterile). Many mating systems have demonstrated 
that when males interact in nature, high-quality males 
usually win the competition for females (Howard et al. 
1997, 1998; Correa and Thiel 2003; Thiel and Correa 
2004). Under some circumstances, however, females may 
not be able to exert their preferences (Hirotani 1994). In 
dung flies and colonial blackbirds, females are only able 
to realize their preferences for high-quality males in low- 
density populations. High-quality males cannot exclude 
other males and low-quality males get access to females 
when population density is high (Borgia 1981; Webster 
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and Robinson 1999; Wong and Candolin 2005). The 
results here shown indicate that high densities of hetero- 
specific males (low quality) interfere with the decision 
process that Drosophila females must make in order to 
choose conspecific over heterospecific males. 

This report demonstrates that at least for these two 
pairs of species of Drosophila, a high relative frequency of 
heterospecific males and a large mating population size 
can affect mating behaviors and lead to increased levels of 
hybridization even in situations in which females have 
access to conspecific males. This might recapitulate the 
situation at the edges of hybrid zones where the relative 
frequency of one of the species is low; in these cases, 
females might be exposed to high ratios of heterospecific 
to conspecific males, which in turn might lead to inter- 
specific matings and increase the likelihood of admixture 
and gene flow. 
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